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Nuclear modification factor Raa{pt) for large transverse momentum pion spectra in Pb + Pb 
collisions at ^/s = 2.76 TeV is predicted within the NLO perturbative QCD parton model. Effect of 
jet quenching is incorporated through medium modified fragmentation functions within the higher- 
twist approach. The jet transport parameter that controls medium modification is proportional to 
the initial parton density and the coefficient is fixed by the RHIC data on suppression of large pt 
hadron spectra. Data on charged hadron multiplicity dNch/drj = 1584 ± 80 in central Pb + Pb 
collisions from the ALICE Experiment at the LHC are used to constrain the initial parton density 
both for determining the jet transport parameter and the 3-l-lD ideal hydrodynamic evolution of 
the bulk matter that is employed for the calculation of RphPhipr) for neutral pions. 
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I. INTRODUCTION 

One important evidence for the formation of strongly 
coupled quark-gluon plasma (QGP) [H-Q in high-energy 
heavy-ion collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) is the observation that the produced dense 
matter is opaque to a propagating parton due to jet 
quenching [5| which suppresses not only the single inclu- 
sive hadron spectra at large transverse momentum 
but also back-to-back high dihadron and 7-hadron 
correlation [9l-[ll|. These observed jet quenching patterns 
in heavy-ion collisions at RHIC can be described by per- 
turbative QCD (pQCD) parton models well ^^M, that 
incorporate parton energy loss [22], [2^ as it propagates 
through dense matter. Since the energy loss or medium 
modification of the effective jet fragmentation functions 
depends on the space-time profile of parton density in 
the medium, any systematic and qualitative extraction 
of the properties of the medium through phenomenolog- 
ical study of jet quenching has to take into account the 
dynamical evolution of the bulk matter (23 - [26| . The ini- 
tial condition for the evolution of bulk matter is either 
provided by model calculation or experimental data on 
the total charged hadron multiplicity. Without any ex- 
perimental data on hadron production in heavy-ion colli- 
sions at the Large Hadron Collider (LHC), all predictions 
for jet quenching 27[ have to rely on theoretical or phe- 
nomenological models on the initial condition for bulk 
matter production and evolution. However, these the- 
oretical and phenomenological predictions for the bulk 
hadron production [23] vary by almost a factor of two 
and; consequently, lead to the same amount of uncertain- 
ties in the predictions for suppression of large transverse 
momentum hadrons in heavy- ion collisions at the LHC. 

Recently, ALICE Experiment at LHC published the 
first experimental data on the charged hadron multi- 
plicity density at mid-ra pid ity in central Pb -\- Pb col- 
lisions at y/s = 2.76 TeV [231. The measured dNch/drj = 
1584 ± 4(stat.) ± 76(sys.) for the top 5% central Pb -t- Pb 



collisions is larger than most of theoretical predictions, 
especially those from the so-called color-glass-condensate 
models [28[ . Such an unexpected large hadron multiplic- 
ity will have important consequences on theoretical pre- 
dictions for jet quenching in Pb+Pb collisions at the LHC 
energies. The predictions for single and dihadron sup- 
pression in the most central Pb+Pb collisions at y/s = 5.5 
TeV by Wang, Wang and Zhang [23| relied on the modi- 
fied HIJING calculation of the charged hadron multi- 
plicity in heavy- ion collisions at the LHC. The modifica- 
tion has now been incorporated into HIJING 2.0 version 
[13 of the original HIJINGl.O model [3l|. The new HL 
JINGH2.0 results [lO,!!! agree well with the first ALICE 
data (28j within experimental errors and theoretical un- 
certainty which is controlled mainly by the uncertainty 
in nuclear shadowing of gluon distribution in nuclei. 

With the first ALICE data [2^ providing more strin- 
gent constraints on the theoretical uncertainty in the bulk 
hadron production, we would like to revisit the predic- 
tions on suppression of single inclusive hadron spectra at 
large px in heavy-ion collisions at LHC. Moreover, we will 
use the 3-l-lD ideal hydrodynamic model for more realis- 
tic description of the bulk matter evolution whose initial 
conditions at LHC are also much better constrained by 
the first ALICE data. 

The rest of the paper is organized as follows. In the 
next section, we will provide a brief overview of the 
pQCD parton model for single inclusive hadron spectra 
and the high-twist (HT) approach to the medium modi- 
fied fragmentation functions. In Sec. HI, we will give a 
brief description of the 3-(-lD ideal hydrodynamic model 
for the bulk evolution with initial condition provided by 
HIJING2.0 model and further constrained by the first 
ALICE data. In Sec. IV, we present predictioirs for the 
nuclear modification factor RAAipr) for neutral pions in 
Pb + Pb colhsions at LHC energy = 2.76 TeV/n and 
discuss about first LHC data for the charged hadrons 
from the ALICE Experiment. We will give a summary 
and discussion in Sec. V. 
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II. PQCD PARTON MODEL AND 
MEDIUM-MODIFIED FRAGMENTATION 
FUNCTIONS 

We will employ the next leading order (NLO) pQCD 
parton model for the initial jet production spectra which 
has been shown to work well for large px hadron produc- 
tion in high energy nucleon-nucleon collisions [33|. We 
use the same factorized form for the inclusive particle 
production cross section in A + B heavy-ion collisions, 
which can be expressed as a convolution of parton distri- 
bution functions inside the nuclei (nucleons) , elementary 
parton-parton scattering cross sections and effective par- 
ton fragmentation functions, 



quenching in dense medium in heavy-ion collisions will 
be described by effective medium modified parton frag- 
mentation functions D^{zh,Q^,E,b,r) within the HT 
approach. 

Within the generalized factorization of twist-four pro- 
cesses, one can calculate the parton energy loss and 
medium modified fragmentation functions of a propagat- 
ing parton in the medium after it is produced via a hard 
process [13, HI] . Within such a high- twist approach, the 
medium modification to the parton fragmentation func- 
tions is caused by multiple scattering inside the medium 
and induced gluon bremsstrahlung. The medium modi- 
fied quark fragmentation functions, 
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where da{ab — >■ cd)/di are elementary parton scatter- 
ing cross sections at leading order (LO) a^. The aver- 
age over the azimuthal angle of the initial jet is implic- 
itly implied in the above equation since we focus on az- 
imuthal integrated single inclusive hadron cross section. 
The NLO contributions include 2 — )• 3 tree level con- 
tributions and 1-loop virtual corrections to 2 — >■ 2 tree 
processes js^] . The nuclear thickness function is normal- 
ized to / d^rt^(r) = A. The nuclear parton distributions 
per nucleon /^/^(xq, r) are assumed to be factorized 
into the parton distributions in a free nucleon fa/N{^i A*^) 
and the nuclear shadowing factor 5*0/^(0;, /^^, r). 



/a/A (a;,M^,r) = S'q/a (a;, /i^, r) 
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where Z is the charge and A the mass number of the 
nucleus. The CTEQ6M parameterization [s^l for par- 
ton distribution functions will be used for nucleon par- 
ton distributions /^^^{x, ^'^). The parton shadowing 
factor Sa/Aix, jJ,'^ ,v) describes the nuclear modification 
of parton distributions per nucleon inside a nucleus 
and can be given by parameterizations [s^. With an 
impact-paramcter-dependent parton shadowing the ef- 
fect of shadowing is the strongest for hard scatterings 
at the center of the transverse plane. But partons from 
these hard scatterings at the center are mostly quenched 
due to parton energy loss and do not contribute to the 
final hadron spectra. Therefore, the effect of parton shad- 
owing is negligible for the final hadron spectra in A + A 
collisions that are dominated by hard scattering close 
to the surface of the overlapping nuclei [s^. We will 
set 5'a/A (a;, /i^, r)=1.0 in the calculation. Effect of jet 
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take a form very similar to the vacuum bremsstrahlung 
corrections that lead to the evolution equations in 
pQCD for fragmentation functions, except that the 
medium modified splitting functions, Ajq^qg(z,£'^) and 
A7,^g,(z,^y) = A75_^gg(l - z,£^) depend on the prop- 
erties of the medium via the jet transport parameter q 

MM, 
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or the average squared transverse momentum broadening 
per unit length, which is also related to the gluon distri- 
bution density of the medium [39l . l40l |. The corresponding 
quark energy loss can be expressed as [1^ [s^ , 
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in terms of the jet transport parameter. Note that we 
include an extra factor ofl — (1 — 2;)/2as compared to 
that used in Refs. [H, |4l| due to corrections beyond the 
helicity amp litude approximation [38| . We refer readers 
to Ref. |26[ for details of the modified fragmentations. 
The fragmentation functions D'^^^^{zc, fi'^) in the vacuum 
arc given by the updated AKK parametrization p^ . 

According to the definition of jet transport parame- 
ter, we can assume it to be proportional to local parton 
density in a QGP and hadron density in a hadronic gas. 
Therefore, in a dynamical evolving medium, one can ex- 
press it in general as 
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where pqgp is the parton (quarks and gluon) density in 
an ideal gas at a given temperature, /(r, r) is the fraction 
of the hadronic phase at any given space and time, 50 
denotes the jet transport parameter at the center of the 
bulk medium in the QGP phase at the initial time tq, 
is the four momentum of the jet and is the four flow 
velocity in the collision frame. We assume the hadronic 
phase of the medium is described as a hadron resonance 
gas, in which the jet transport parameter is approximated 



as. 
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where pM and pB are the meson and baryon density in 
the hadronic resonance gas at a given temperature, re- 
spectively, Pn ^ no Ki 0.17 fm^'^ is the nucleon density in 
the center of a large nucleus and the factor 2/3 accounts 
for the ratio of constituent quark numbers in mesons and 
baryons. The jet transport parameter at the center of a 
large nucleus qm has been studied in deeply inelastic scat- 
tering (DIS) |43l - l4a | . We use a recently extracted value 
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QN « 0.02 GeVVfm from the HERMES [46| experi- 
mental data. The hadron density at a given temperature 
T and zero chemical potential is 
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where r;^ = ± for meson (M) /baryon (B). In the paper, 
we will include all hadron resonances with mass below 1 
GeV. 



III. 3+lD IDEAL HYDRODYNAMIC 
EVOLUTION OF BULK MATTER 

In the model for medium modified fragmentation func- 
tions as described in the last section, one needs informa- 
tion on the space-time evolution of the local tempera- 
ture and flow velocity in the bulk medium along the jet 
propagation path. We will use a full three-dimensional 
3-l-lD ideal hydrodynamics [Ull^ in our calculation to 
describe the space-time evolution of the bulk matter in 
heavy-ion collisions. 

We solve equations of energy-momentum conservation 
in full 3-l-lD space (r, x, y, 77.,) under the assumption that 
local thermal equilibrium is reached at an initial time 
To =0.6 fm/c and maintained thereafter until freeze-out. 
Here r, rjg, x and y are proper time, space-time rapid- 
ity, and two transverse coordinates perpendicular to the 
beam axis, respectively. Ideal hydrodynamics is charac- 
terized by the energy-momentum tensor. 



(9) 



where e, P and arc energy density, pressure and lo- 
cal four velocity, respectively. We neglect the finite net- 
baryon density which is supposed to be small both at 



RHIC and LHC energies. For the quark-gluon plasma 
(QGP) phase at high temperature [T > = 170 MeV), 
we use the equation of state (EOS) of a relativistic mass- 
less parton gas (u, d, s quarks and gluons) with a bag 
pressure B: 



P=\{e-AB) 



(10) 



The bag constant is tuned to be = 247 McV to match 
the pressure of the QGP phase to that of a hadron res- 
onance gas at critical temperature Tc = 170 MeV. Be- 
low the critical temperature T < Tc, a hadron resonance 
gas model including all hadrons up to the mass of the 
A(1232) is employed for the EOS . The hadron reso- 
nance gas EOS employed in this study implements chem- 
ical freeze-out at Teh = 170 MeV [43 so that evolution 
of chemically frozen, but thermally equilibrated hadronic 
matter is described. In the calculation of parton energy 
loss and medium modified fragmentation functions, jets 
propagate along a straight path in the evolving medium 
until the decoupling of the medium at r^oc — ^00 MeV. 

For the initial condition of the longitudinal flow ve- 
locity, Bjorken's scaling solution [i^ is employed. The 
initial entropy distribution in the transverse plane is pro- 
portional to a linear combination of the number density of 
participants, Ppart, and that of binary collisions, pcoii |50j . 
The proportionality constant and the fraction of soft and 
hard components are so chosen that the centrality depen- 
dence of charged particle multiplicity agrees with the ex- 
perimental data at RHIC ^ and HIJING2.0 results at 
the LHC energy which agree with the recent ALICE 
data for the most central Pb + Pb collisions at y/s = 2.76 
TeV [15. Shown in Fig. [1] is the hydrodynamic calcu- 
lation (solid square) of the charged hadron multiplicity 
density at mid-rapidity as compared with the HIJING2.0 
result (solid circle). 

We will use these hydrodynamic solutions to provide 
values of temperature and flow velocity along the tra- 
jectory of parton propagation for the evaluation of local 
jet transport parameter in Eq. ([6]) for the calculation of 
medium modified fragmentation functions. The hydro- 
dynamic results obtained in this way such as tempera- 
ture and energy density at each space-time position are 
publicly available [5l|. Note that these hydrodynamic 
results provide the evolution of the space-time profiles 
of jet transport parameter and gluon density relative to 
their values at the center of overlapped region of dense 
matter (r = 0) in the most central collisions (& = 0) at 
an initial time tq . The normalization of the initial values 
at To and r = 0,6 = is included in the value of go- 
It should be proportional to the initial parton density, 
which in turn is determined by fitting the final charged 
hadron multiplicity density in mid-rapidity (see Fig. [T]) 
of the hydrodynamics results. Therefore, approximately. 
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We determine the coefficient of the above relation by fit- 
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FIG. 1: (Color online) Charged particle multiplicity density at 
midrapidity per participant pair as a function of the number 
of participants in Pb + Pb collisions at ^/snn = 2.76 GeV 
from HIJING 2.0 (solid circles) [s^l and ideal hydrodynamics 
(solid square). 



ting the experimental data at the RHIC energy for the 
most central Au + Au collisions. Then the energy de- 
pendence and the impact-parameter dependence will be 
completely determined by the measurement or calculated 
values of charged hadron multiplicity density dNch/dr]. 



A. Numerical results 

We use a NLO Monte Carlo based program [11] to cal- 
culate the single hadron spectra in our study. In this 
NLO program, the factorization scale and the renormal- 
ization scale are chosen to be the same (denoted as /i) 
and are all proportional to the transverse momentum of 
the final hadron p^- We choose /i = 1.2pT with which the 
calculated single inclusive pion spectra for p + P collisions 



in NLO pQCD agree weh with RHIC data [171. We will 
use the same scale in A + A collisions at LHC energies. 

With space-time profile of the gluon density provided 
by the hydrodynamic evolution equations we calculate 
medium modified fragmentation functions which are used 
then to calculate the suppression factor (or nuclear mod- 
ification factor) for large px hadron spectra in heavy-ion 
collisions [52l. 
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where N^,^{b) = J d'^rtA{r)tB{\b - r\). The fixed value 
of impact-parameters in the calculation of the spectra 
and the modification factor are determined through the 
Glauber geometric fractional cross sections for given cen- 
trality of the heavy-ion collisions. 

In the HT approach, we have determined QoTq ~ 
0.54 — 0.63 GeV^ from the experimental data on pion 



spectra in the most 0-10% central Au + Au collisions at 
•y/s = 0.2 TeV [H]. We assume that the jet transport 
parameter is proportional to the initial parton density 
or the transverse density of charged hadron multiplicity 
in mid-rapidity [Eq. [TT] . With the new ALICE data on 
charged particle pseudo-rapidity density at mid-rapidity 
dNch/dri = 1584 ± 4(stai.) ± 7&{sys.) [H] in the most 
central 5% Ph + Pb colhsions at a/s = 2.76 TeV ver- 
sus dNch/drj = 687 ± 37 for 0-5% Au + Au collisions 
at \/s = 0.2 TeV [sl], we obtain the extrapolated value 
goTo = 1.0- 1.4 GeV^ for Pb + Pb collisions at = 2.76 
TeV. Shown in Fig. [2] is the predicted nuclear modifica- 
tion factor for pion spectra in the 0-5% central Pb + Pb 
collisions at = 2.76 TeV. The suppression factor in- 
creases with pt partially because of the energy depen- 
dence of parton energy loss and partially because of the 
less steep initial jet production spectra [H^. This trend 
is similar to almost all LHC predictions by many other 
parton energy loss model calculations H^] ■ 

We also show in Fig. [5] the recently published AL- 
ICE data [m (filled square) on the suppression factor 
for charged hadrons in the most 0-5% central Pb + Pb 
collisions at the LHC energy y/s = 2.76 TeV. Since there 
are no experimental data on charged hadron spectra in 
p+p collisions at the LHC energy = 2.76 TeV, ALICE 
data on the suppression factor were obtained with p+p 
spectra interpolated from experimental data at Tevatron 
energies. The histograms in Fig. [2] represent the errors on 
the suppression factor from the uncertainty in the inter- 
polation. Charged hadrons contain significant fraction of 
protons and anti-protons which could have nonnegligible 
contributions from parton recombination |56l - l58l |. Be- 
cause of the abundance of jet production at the LHC en- 
ergy, recombination among these hard partons becomes 
possible and therefore contribute to hadron, especially 
baryon spectra at high px- One therefore should take 
into account the contribution from hard parton recombi- 
nation in the calculation of the final charged hadron spec- 
tra, which could push the suppression factor for charged 
hadrons higher than that for pions at large pr ■ 

Because of the increased initial parton density which 
we assume to be proportional to the final hadron multi- 
plicity density, the initial jet transport parameter go in 
Pb + Pb colhsions at the LHC energy =2.76 TeV are 
more than twice larger than that in Au + Au collisions at 
the RHIC energy y/s = 0.2 TeV. The hadron suppression 
factors in Pb + Pb collisions at LHC at moderate trans- 
verse momentum px < 20 GeV/c are therefore about 
twice smaller than that at RHIC [26|. Shown in Fig. El 
arc the scaled jet transport parameter g(r, r)/q(0, tq) as 
a function of t — tq which is related to the parton and 
hadron density [Eq. ([6])] as given by the hydrodynamic 
evolution. The kink in the time dependence is caused by 
the first-order phase transition assumed in the hydrody- 
namics evolution as the EOS of the dense matter. For 
most part of the evolution history, the scaled jet trans- 
port parameters are very similar at RHIC and LHC en- 
ergies. Therefore, the increased hadron suppression at 
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I 0-5% Pb+Pb^;i° \ls=2.76TeV q^T„=1.0-1.405GeV^: 




FIG. 2: (Color Online) Nuclear modification factor at mid- 
rapidity for neutral pion spectra in the most — 5% central 
Pb + Pb collisions at ^ = 2.76 TeV, using the HT modified 
fragmentation functions with qoro = 1.0—1.4 GeV^. The data 
points (filled square) are for charged hadrons in the same cen- 
tral Pb + Pb collisions from ALICE Experiment [5^1 with the 
histogram representing systematic errors due to uncertainty 
in the p + p spectra at = 2.76 TeV from interpolation. 



predict the suppression factor for the hadron spectra in 
Pb+Pb coUisions at = 2.76 TeV within the HT model 
for medium modified fragmentation functions. Because 
of the increased initial parton density of about a factor of 
two, and the longer life-time of the dense matter or the 
duration of jet quenching, the hadron spectra are found 
to be suppressed more at LHC than at RHIC. However, 
because the energy dependence of the parton energy loss 
and the less steep initial jet spectra, the suppression fac- 
tors will increase with px more strongly at LHC than at 
RHIC. 

Because of the increased number of jet production in 
heavy-ion collisions at LHC energies, there are increased 
possibility of larger pT hadron production from recom- 
bination of parton showers from independent jets. This 
production mechanism will be more important than the 
shower-thermal and thermal-thermal parton recombina- 
tion that have been considered more relevant in heavy-ion 
collisions at the RHIC energy [56, ,58J . Such contributions 
from jet-jet parton recombination will likely increase the 
hadron yield at moderate pt and increase the values of 
the suppression factor Raa at the LHC energies. 



LHC is caused mainly by the overall increase of the ini- 
tial parton density. The increased initial parton density, 
however, will also increase the life-time of the dense mat- 
ter through-out the phase transition and hadronic phase. 
This will also contribute to the increased suppression of 
hadron spectra at the LHC as compared to at RHIC. 

IV. CONCLUSIONS 

We used the new ALICE data on charged hadron mul- 
tiplicity density at mid-rapidity in central Pb -\- Pb colli- 
sions at the LHC energy ^/s = 2.76 TeV [2^ to estimate 
the initial jet quenching parameters in Pb + Pb colli- 
sions at the LHC and the initial condition for the hydro- 
dynamic evolution of the bulk matter. With the initial 
values of the jet transport parameter and the initial con- 
dition for hydrodynamic evolution of the bulk matter, we 
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